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EXTERNAL SHEAR IN QUADRUPLY IMAGED LENS SYSTEMS 
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ABSTRACT 

We use publicly available N-body simulations and semi-analytic models of galaxy formation to estimate the 
levels of external shear due to structure near the lens in gravitational lens systems. We also describe two selection 
effects, specific to four-image systems, that enhance the probability of observing systems to have higher external 
shear Ignoring additional contributions from "cosmic shear" and assuming that lens galaxies are not significantly 
flattened, we find that the mean shear at the position of a quadruple lens galaxy is 0.11, the rms shear is ^0.15, 
and there is a ^45% likelihood of external shear greater than 0. 1 . This is much larger than previous estimates and 
in good agreement with typical measured external shear. The higher shear primarily stems from the tendency of 
early-type galaxies, which are the majority of lenses, to reside in overdense regions. 



1. INTRODUCTION 

With the benefit of hindsight it can be argued that even the 
second gravitational lens discovered, PGl 115-1-080 (Weymann 
et al. 1980), foreshadowed what has turned into an embar- 
rassment of riches - a superabundance of quadruply imaged 
quasars. But only with the advent of systematic lens surveys 
(King and Browne 1996; Rusin and Tegmark 2001) has it be- 
come clear that the high ratio of quadruple to double systems is 
not an artifact of observational selection and therefore presents 
a genuine challenge to our understanding of lensing of galax- 
ies. Furthermore, individual fits on a system-by-system basis 
often require large amplitude (0.1-0.3) external shear (Hogg 
and Blandford 1994; Schechter et al. 1997; Kneib, Cohen, 
and Hjorth 2000; Fischer, Schade, and Barrientos 1998), sig- 
nificantly higher than values of 0.02-0.05 expected (Keeton, 
Kochanek and Seljak 1997; henceforth KKS) from large scale 
structure or nearby galaxies. 

There are three factors that will probably have some part to 
play in the ultimate resolution of these problems: galaxy ellip- 
ticities, shear due to random superpositions of mass along the 
line of sight, and shear due to structures that are associated with 
the lens galaxy. In this paper we concentrate on the shear from 
associated structures. 

The relative importance of tides and ellipticity has been con- 
sidered by KKS. In computing the expected tidal shear they 
consider three contributions: a) random shear due to unassoci- 
ated foreground and background structure, b) the effect of asso- 
ciated galaxies, through the two point correlation function, and 
c) the effect of clusters of galaxies, with a term proportional to 
the number density of clusters. 

In the present paper we take a different tack to estimate the 
expected tidal shear due to nearby structures. We use the GIF 
Project recipe (Kauffmann et al. 1999) for galaxy formation 
within a cold dark matter (CDM) simulation to compute the 
shear expected along random lines of sight, in the directions of 
galaxies, and specifically in the directions of early-type galax- 
ies, which appear to be the predominant type of lens galaxy. 
In this way, the effects of correlated structure are naturally in- 
cluded, without artificial distinctions between different types of 
structures, and we can also include the clustering properties of 
different types of galaxies. This is similar to recent work by 

' http: //www.MPA-Garching.MPG.DE/GIF 



White, Hernquist, and Springel (2001), where a high resolu- 
tion hydrodynamical simulation was used to find typical values 
of external shear at the positions of typical galaxies, rather than 
specifically at the positions of massive early-type galaxies. 

In § 2 we outline our methods for using the GIF simulations 
to estimate the effects of correlated structure on the statistics of 
external shear in gravitational lens systems. In § 3 we describe 
two selection effects which will further increase the typical ex- 
ternal shear measured in quadruple gravitational lens systems; 
high shear systems have a larger cross-section for quadruple 
lensing (which is partly offset by "magnification bias") and re- 
gions of high external shear are likely to be regions of higher 
convergence, due to large scale structure. We finish with a 
discussion of the consequences of our calculations for various 
problems associated with quadruply imaged systems. 

2. LARGE SCALE STRUCTURE AND TYPICAL SHEAR 

Studies of cosmic shear have noted rms shear of a few per- 
cent (Bacon et al. 2001 ; Van Waerbeke et al. 2001 ; Wittman 
et al. 2000) from studies of weak distortions of background 
galaxies. This result can not be directly applied to estimates of 
typical shear values for lens systems because a strong gravita- 
tional lens is not at a random position in the sky but is instead at 
the position of a fairly massive galaxy, typically an early-type 
(elliptical or SO) galaxy (Keeton, Kochanek, and Falco 1998). 
Such galaxies are known to be located preferentially in over- 
dense regions, where it would also be expected that the typical 
shear would be higher. 

To estimate the effects of large scale structure that is cor- 
related with the lens galaxies, we used the publicly available 
simulations from the GIF Project'. The simulations provide the 
dark matter distribution as well as estimates of the positions, ve- 
locities, luminosities and colors of galaxies (Kauffmann et al. 
1999). The galaxy information is derived from a semi-analytic 
model of galaxy formation. The simulation box was 141.3 
/i"'Mpc (comoving) on a side, the 256^ particles each had a 
mass of 1.4x IO^'/i'^Mq, the gravitational softening length was 
20/!~'kpc, and the cosmological model was a flat universe with 
f7,„ = 0.3, (78 = 0.9, /i = 0.7, r = 0.21. 

We made projected mass maps of the outputs at z = 0.42 at a 
resolution of 50 /i"' kpc and used these maps to generate shear 
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maps, assuming a source at = 3. The critical surface den- 
sity for lensing for our assumed cosmology is approximately 
2800/iMqPC"^. At this redshift the box size corresponds to 
Az ^ 0.05 and our comoving 50 /;"' kpc resolution element cor- 
responds to an approximate angular size of 8", much larger 
than the typical Einstein radius of a galaxy 1")- 

The relation between convergence (mass surface density in 
units of the critical surface density) and shear are algebraic in 
Fourier space (e.g., Kaiser and Squires 1993), so we used an 
FFT method to generate shear maps from convergence maps. 
The Fourier transforms of the convergence map, k, and the two 
shear components, 7,, are simply related: 



7i(m,v)= 



k(u,v) ; 72(m,v): 



-2uv 



k(u,v) , (1) 



where a tilde denotes a quantity in the Fourier domain and u and 
y are Fourier space variables (i.e., wavenumbers). After mak- 
ing the convergence map from the particle data of the N-body 
simulation, the map was Fourier transformed, the algebraic op- 
erations were performed to construct the shear components in 
the Fourier domain and the Fourier shear maps were inverse 
transformed to produce shear image maps. 

With the GIF catalog, we measured the shear at all positions 
that were marked as having massive early-type galaxies. Our 
criteria for being a massive early-type galaxy were B-V > 0.9 
and Msiar > 10"/;~'Mq. We also measured the shear at the 
position of all galaxies in the GIF catalog, as well as at a com- 
parable number of random points in the map. 

The measured shear at random points, shown in figure 1, was 
typically less than 1 %, as would be expected from a box with 
redshift extent of approximately Az ^ 0.05. The few percent 
cosmic shear is built up from having many volumes of this ex- 
tent along the line of sight. Typical shear values around galaxies 
were noticeably higher, while the mean shear at the position of 
bright early-type galaxies was yet 
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Fig. 1 . — Effects of galaxy bias on typical sheai' in lens systems. Black solid 
curves show shear distributions at position of early-type galaxies, red dashed 
curves show shear distributions at position of all galaxies, blue dotted curves 
show shear distributions at random positions in image. Three curves for each 
type are for three projections of the simulation volume (along the simulation's 
x,y and z axes). 

higher, 0.058, with an rms shear of 0.071 and an extensive 
tail past 7 ~ 0.1. The probability of shear greater than 0.1 is 



0.16. This is a consequence of early-type galaxies being pref- 
erentially found in or near clusters of galaxies. Previous esti- 
mates (e.g., KKS) took into account that galaxies were corre- 
lated with other galaxies in estimating external tidal effects, but 
ignored the preferential location of lenses inside larger struc- 
tures in estimating the contribution from large scale structure. 
This is an important consideration for two reasons. More mas- 
sive galaxies, which have a higher cross-section for producing 
gravitational lens systems, are strongly correlated with large 
scale structure. Furthermore, early-type galaxies are even more 
highly correlated with large scale structure than late-type galax- 
ies. 

Thus we see that the material correlated with the lens galax- 
ies can contribute a significant amount of shear. The results pre- 
sented here are probably an underestimate of the likelihood of 
high shear events, as the finite box size results in missing long- 
wavelength power. Other material along the line of sight will be 
much less correlated, but will add a significant amount of rms 
shear (KKS find typically ~ 3% for this cosmological model). 
The typical shear at the position of a massive early-type galaxy 
is therefore likely larger than ^1%. Note that this is applicable 
to all gravitational lens systems (and early-type galaxies in gen- 
eral) and that we have not introduced any distinctions between 
two-image and four-image systems. 

In making the maps smoothing was done on a scale of ap- 
proximately 50 /;"' kpc, reducing the effects of nearby galax- 
ies. There is a trend in our results for increasing shear with 
increasing resolution, suggesting a non-negligible amount of 
shear from near neighbors, as also suggested by White, Hern- 
quist, and Springel (2001); with a smoothing length of 200 
kpc the mean shear at the position of early-type galaxies 
drops to 0.03. The larger smoothing primarily reduces the prob- 
ability of high shear values, perhaps suggesting that higher res- 
olution would lead to an even larger probability of high shear. 
In the case of the gravitational lens system B 1422-1-231 (Hogg 
and Blandford 1994) it appears that most of the shear (> 0.2 
in this case) does in fact come from neighbors closer than 50 
kpc. 

When all particles within a spherical radius of 50 /;"' kpc 
(comoving) from the centers of early-type galaxies were re- 
moved before creating the initial mass map the statistics were 
nearly unchanged; for example the mean shear decreased by 
only 0.003. When this radius of excision was doubled the mean 
shear was reduced to 0.045, a fractional reduction of more than 
20%. This agrees with both the above discussion on smoothing 
and the work of White, Hernquist, and Springel (2001), and 
points to the importance of substructure that is associated with 
the lens galaxy but in the outer regions of the dark matter halo. 

3. FURTHER SELECTION EFFECTS 

The previous section demonstrated that not all positions on 
the sky are equally likely to host a bright early-type galaxy; 
it is also true that not all positions on the sky hosting bright 
early-type galaxies are equally likely to host a quadruple grav- 
itational lens system. The cross-section for quadruple lensing 
by a galaxy can be enhanced by external shear, while the mag- 
nification statistics are affected by both external shear and con- 
vergence from large scale structure. A bright early-type galaxy 
in a region of high shear and high convergence has a higher 
probability of being observed to be a quadruple gravitational 
lens. 

In this section, we calculate the relative number of quadruple 
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images expected in samples, taking into account three effects: 
1) cross-section for quadruple lensing increases with increas- 
ing shear; 2) magnification for quadruple systems decreases 
with increasing shear; 3) magnification increases with increas- 
ing convergence. 

The cross-section for producing a quadruple lens is simply 
the solid angle subtended by the part of the source plane that is 
able to produce four observable images (and a fifth highly de- 
magnified, thus very difficult to observe, image) for a given lens 
position; i.e., the cross-section is the area inside the "diamond 
caustic." For a singular isothermal sphere (SIS) lens model in 
the presence of external shear, the cross-section for producing 
a quadruply imaged system scales as 7^/(1 - 7^)(Finch et al. 
2002). The probability of an object producing four images 
therefore scales roughly as 7^, increasing the probability of 
galaxies in the presence of high external shear being in a catalog 
of observed lenses. 

The magnification is important through the effects of "mag- 
nification bias" (Turner, Ostriker, and Gott 1984). The number 
of quasars increases toward lower unmagnified flux; high mag- 
nification increases the possible number of observable sources. 
The effects of magnification bias are often in the opposite direc- 
tion of the change in cross-section. Calculation of the expected 
number of lens systems requires the number of sources per unit 
sky area (in the source plane) that could be sufficiently magni- 
fied to enter the sample. 

The number of sources per source plane solid angle expected 
above a given observed flux limit 5m,„, in a region of magnifica- 
tion /i, and assuming a power-law unmagnified flux distribution 
dN/dS = NoiS/Sor" is 

N(>S,„i„)= / Noi — )dS. (2) 

Employing a change of variable, this reduces to simply 
times the number per source plane solid angle expected in the 
absence of magnification. For a single lens system, there is a 
distribution of magnifications in the source plane, so the fac- 
tor by which magnification increases the number of observ- 
able quasars is then < >, where the average is over the 
source-plane distribution of magnifications. For the special 
case of a = 2, close to the observed value (Rusin and Tegmark 
2001), the factor by which the number of observable sources 
per source plane solid angle is increased (often called "magni- 
fication bias") is equal to the mean magnification. 

The mean magnification of a quadruple system for an SIS 
model in the presence of external shear (Finch et al. 2002) 
scales as 7~'(1 ~7^)~'- Magnification bias, in the presence of 
external shear, therefore has the effect of reducing the number 
of observed lens systems, but by a factor that is smaller than that 
by which the cross-section is increased. The net effect of exter- 
nal shear is to enhance the number of quadruple lenses, and we 
would expect instances of high external shear to be overrepre- 
sented in catalogs of quadruple gravitational lenses relative to 
massive early-type galaxies in general or two-image lenses. 

The effect of additional convergence k,, is to scale the mag- 
nifications by a constant factor (1 - k^)"^ and to enhance the ef- 
fects of external shear, causing an "effective shear" (Paczynski 
1986) Jeff = 7/(1 ~ Ks). In practice, for the shear and conver- 
gences due to large scale structure the difference between 7 and 
Jeff is negligible. 

If we approximate the effects of magnification bias as scaling 
the fraction of lenses by the mean magnification (i.e., assume 
a = 2), the net effect (including cross-sections and magnifica- 



tion) of external shear and additional convergence is that the 
probability of a system being observed (assuming a SIS lens 
with external shear) to be a quadruple lens is proportional to 
7f//(l - KsT^- The tendency of early-type galaxies 

to reside in overdense regions, where both external shear and 
additional convergence are likely, leads to a significant selec- 
tion effect in the direction of enhancing the typical shear at the 
position of a quadruple image lens system. 

This is shown in figure 2, where we show the effects of this 
selection on the shear (true shear, not effective shear) expected 
for a typical four-image lens system, assuming that the lens can 
be modeled as an SIS with external shear. As the input shear 
distribution we used the results from the previous section, us- 
ing only the contribution from the material in the simulation 
volume and not including the contribution from other mate- 
rial along the line of sight. The shear distribution was then 
constructed, with each galaxy position weighted by the appro- 
priate quadruple lens probability. The mean convergence at 
each galaxy position was measured using the original mass map 
smoothed on a scale of 300 /i"' kpc to ensure that the additional 
convergence was truly from large scale surrounding material. 
This will underestimate the extra convergence expected at the 
position of lens galaxies, but significantly reduces shot noise as 
well as ambiguity concerning what should be considered part 
of the lens galaxy. 

Including these selection effects, the mean shear expected 
from only the immediately surrounding material is 0. 1 1 and the 
rms shear is 0.15, with a ~45% probability of 7 > 0.1. There- 
fore, it is to be expected that measured shear in typical four- 
image systems should be > 0.1. 
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Fig. 2. — Effects of observational biases on observed average shear in four- 
image lens systems, assuming a lens model with a SIS and external shear. Blue 
dotted curve is shear distribution at position of early-type galaxies (same as the 
solid curves in figure 1), red dashed curve shows selection effects due to ex- 
ternal shear only, black solid curve shows combined effects of external shear 
and extra convergence. Each curve shows the average histogram of the three 
viewing projections. 



We have assumed that galaxy ellipticity is not important, 
but ellipticity of the lens galaxy can significantly reduce the 
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enhancement of the cross-section due to external shear. For 
galaxies more elliptical than an axis ratio of about 0.7:1 the 
cross-section is nearly unaffected by external shear (Rusin 
and Tegmark 2001). Observed samples of early- type galaxies 
in nearby clusters Jorgensen, Franx, and Kjaergaard (1995) 
contain many significantly flattened objects, so if typical 
lens galaxies are morphologically similar to typical early-type 
galaxies (of all masses) in galaxy clusters, then the selection 
effects outlined in this section wiU be reduced. 

4. OUTSTANDING ISSUES IN QUADRUPLY IMAGED SYSTEMS 

Previous studies have concluded that external shear is not 
a viable solution to the apparent excess of quadruples (Rusin 
and Tegmark 2001 and references therein) but on a system-by- 
system basis it is apparent that external shear is significantly 
higher in lens systems than had been predicted. The calcu- 
lations in the previous sections demonstrate that the external 
shear expected in a typical quadruple lens system is roughly 
0.10, with a large tail extending to ~ 0.3, much larger than pre- 
viously predicted. The effects of large scale structure uncorre- 
lated with the lens, structures within 50 kpc of the lens, and 
substructures within the lens galaxy wiU increase typical shear 
values further. 

Using the compilation of lens models from the CASTLES 
web site^, the mean reported external shear (for all lens sys- 
tems, not just quadruple systems) is 0.10, with an rms of 0.12. 
This should be treated as anecdotal evidence for typical shear 
of ~0.1; several models have a non-positive number of degrees 
of freedom, the reported uncertainties have been ignored in 
our calculations (i.e., not inverse variance weighted) and many 
of the fits have been performed holding the galaxy ellipticity 
fixed. We have shown that such a distribution of shears can 
comfortably be accommodated within the current understand- 
ing of galaxies and their environments. The observation that 
lens galaxies are predominantly early-type galaxies is fully con- 
sistent with the observed external shears of ~ 0.1. As would be 
expected, the observed shear in quadruple-image systems ap- 
pears to be systematically larger than for double-image lenses. 

To solve the excess of four-image systems compared to two- 
image systems, typical external shear on the order of 0.3 is re- 
quired (KKS; Finch et al. 2002). We are in agreement with 
previous work that such high shears are not expected from cur- 
rently viable cosmological models. A combination of external 
shear and galaxy ellipticity is not efficient for increasing the 
fraction of quads, as they are most likely not correlated and 
therefore do not add coherently. The main effect of external 
shear in the presence of a distribution of ellipticities is to in- 
crease the fraction of quadruple images only in the low elliptic- 
ity regime. 

Other possible solutions to the problem of too many quadru- 
ple images can be enhanced by the addition of external shear. 
For example, adding a core to the mass profile or making the 
central density profile less steep can enhance the fraction of 
quadruple images (Rusin and Tegmark 2001). The additional 
convergence due to large scale structure will introduce a moder- 



ate selection effect, as the image separations will be increased, 
making it more likely that small splittings could be resolved. 
However, it does not appear that external shear can by itself 
solve the problem of the apparent excess of four-image systems. 

5. DISCUSSION 

We have shown that the relatively high amounts of exter- 
nal shear found in most lens systems naturally arise from non- 
linear large scale structure in the vicinity of the lenses. Early- 
type galaxies are systematically in more overdense regions and 
therefore would be expected to experience higher amplitude 
tidal fields. 

The fact that quadruple lens galaxies are mainly early-type 
galaxies is enough to explain the observed values of external 
shear, previously considered high. The dominance of early-type 
galaxies in lens samples arises simply from the bias in lens sam- 
ples toward more massive galaxies, which are more likely to be 
early-type galaxies (Fukugita and Turner 1991; Maoz and Rix 
1993; Kochanek 1993). 

This amount of shear expected from nearby structures is not 
sufficient to, by itself, solve the "quad problem" (the apparent 
relative super-abundance of four-image systems) but the effect 
goes in the right direction and should add with other effects, 
such as making density profiles in the centers of galaxies less 
steep. Paying careful attention to the preferred positions of 
galaxies relative to typical places in the universe, as we have 
done here, will almost certainly have an important role in un- 
derstanding these observations. 

Typical shear values are sufficiently high that external shear 
should be of comparable or greater importance in lens mod- 
eling to internal shear (galaxy ellipticity), in agreement with 
model fits to known lenses. We do not expect the introduction 
of galaxy ellipticities to significantly change our conclusions, 
but it would be a natural direction for future work. 

Selection effects are important in gravitational lens systems. 
The results presented here suggest that our knowledge of struc- 
ture formation and galaxy environments can be usefully applied 
to understanding the observed statistics and characteristics of 
lens systems. As numerical simulations of large scale structure, 
models of galaxy formation, and gravitational lens catalogs all 
continue to improve we can expect the importance of an inte- 
grated approach to only increase. 
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